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ABSTRACT 


A  new  finite  element  modeling  is  presented  to 
investigate  the  static  and  dynamic  behavior  of  laminated 
composite  beams  with  partial  delamination.  In  this  study,  a 
newly  developed  rectangular  be'?)*  element  is  used.  The  element 
has  lateral  and  axial  displacen  ^nt.  as  degrees  of  freedom  but 
no  rotation.  For  simplicity,  linear  shape  functions  are  used 
for  the  beam  element.  As  a  result,  the  element  has  six  degrees 
of  freedom,  four  of  which  are  the  axial  displacements  at  the 
corner  points  and  two  are  the  lateral  displacements  at  the 
ends.  In  addition,  contact -impact  conditions  are  applied  to 
the  finite  element  modeling  to  avoid  overlapping  of  the  upper 
auid  lower  portions  of  a  delaminated  section.  The  numerical 
study  shows  that  depending  on  existence  of  an  embedded 
delamination  crack  and  its  size,  the  response  is  different  for 
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I .  INTRODUCTION 


Laminated  composite  structures  provide  improved  mechanical 
properties  such  as  high  strength,  stiffness,  toughness  and 
high  temperature  performance  [Ref.  1]  .  However,  these 
structures  have  more  complicated  failure  mechanisms  and  modes 
compared  to  those  of  conventional  metallic  structures.  One  of 
the  common  failure  mode  is  inter-delamination. 

The  objectives  of  this  study  are  to  develop  a  finite 
element  model  for  an  analysis  of  laminated  composite  beam 
structures  containing  an  inter-delcunination  crack  and  to  find 
a  methodology  to  detect  an  embedded  interlaminar  delamination. 
For  delaminated  beams,  a  short  impulse  load  is  applied  to  find 
out  whether  the  dynamic  response  varies  significantly  compared 
to  the  dyneunic  response  of  structures  without  a  crack. 

In  this  study,  a  newly  developed  finite  element  model 
[Ref.  2]  is  used  for  the  simulation  of  beams  containing  an 
inter-delamination  crack.  A  recteuigular  element  that  has  six 
degrees  of  freedom  is  used.  The  element  allows  axial 
displacements  at  its  four  comer  points  and  lateral 
displacements  at  its  two  ends.  Using  this  element,  static  and 
dynamic  analyses  of  laminated  composite  beams  containing  a 
local  delamination  crack  are  conducted.  First  of  all,  static 
deflections  as  well  as  natural  frequencies  and  mode  shapes  of 
beams  without  delamination  are  calculated  to  check  the 
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II.  FINITE  ELEMENT  FORMULATION 


A.  FINITE  ELEMENT  MODELING 

A  newly  developed  rectangular  finite  element.  Figure  1,  is 
used  to  investigate  the  static  and  dynamic  behavior  of 
laminated  composite  beams.  The  element  has  six  degrees  of 
freedom  and  allows  axial  cind  lateral  displacements  but  no 
rotation.  There  are  axial  displacements  at  the  four  corner 
points  and  lateral  displacements  at  the  two  ends  of  the 
element . 


Figure  1  Four  Noded  Rectangular  Element  with  Six  Degrees  of  Freedom. 


In  Figure  1  euid  the  subsequent  formulation,  u  represents 


the  axial  displacements  at  the  corner  points  and  v  represents 
the  lateral  displacements  at  the  ends.  The  subscripts  “1"  and 
"2"  refer  to  the  left  and  right  ends  while  the  superscripts 
"t“  and  "b"  indicate  the  top  and  bottom  sides  of  the  element, 
respectively . 

The  displacement  model  of  the  element  is 

(1) 

where  [N]  is  the  matrix  of  shape  functions  and  {d*}  is  a 
vector  of  nodal  displacements .  The  axial  displacement  is 
assiimed  to  vary  linearly  along  both  axial  and  lateral 
directions .  It  can  be  written  as 


u  (x,  y)  (X)  [ja^  (y)  u/+i^  (y)  u/] 

(X) (y)  (x) (y)  ^ ^  Uj" 

The  lateral  displacement,  which  is  assumed  to  be  constant 
through  the  thickness  of  the  element,  varies  linearly  along 
the  cixial  direction  and  can  be  written  as 


v(x)  =  Z)w^(x)yj 

i-l 

=W  (X)  Vj+Jjr  (X)  Vj 


(3) 


Here,  Ni  and  Hi  are  the  linear  interpolation  or  shape  functions 
in  the  axial  and  lateral  directions.  The  beam  element  may  use 
a  higher  order  shape  function  for  the  axial  direction,  i.e. 
Ni,  and  the  linear  shape  function  for  the  lateral  direction, 
i.e.  Hi.  However,  in  this  study,  the  linear  shape  fxinction  is 


used  for  both  Ni  and  Hi  for  siit^Jlicity .  That  is. 


Wi(x)=l-- 
^  a 


(4) 


where  a  and  b  are  the  length  and  height  of  che  beam  element, 
respectively . 

For  simplicity,  notations  N^,  Nj,  Hi  aind  H2  will  be  used 
instead  of  Ni  (x)  ,N2  (x)  ,Hi  (y)  and  H2(y)  in  the  following 
derivation. 

For  two  dimensional  elasticity  problems  the  strain  vector 
can  be  expressed  as  [Ref.  3] 

-  h' 


d/dx  0 
0  d/dy\ 
d/dy  a/sJ 


(5) 


where  {d*}  is  the  vector  of  nodal  displacements  and  [B]  is  the 
matrix  that  relates  strains  to  the  nodal  displacement  vector . 
Nozrmal  strains  can  be  written  as 


*  dx  dx  dx 

£  =-^=0 
dy  ^ 


(6) 


and  the  shear  strain  is 
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dy 


du  ^  dv 
dx 


dH^ 

~dy 


NX 


dH,  .  dlL  .  dN  dN 


3x^2 


(7) 


The  element  stiffness  matrix  can  be  obtained  by  minimizing 
the  total  strain  energy  which  contains  both  bending  and  the 
transverse  shear  energry.  This  minimization  yields  the 
following  element  stiffness  matrix 

lk*}  =  [k^n*[k^n  (8) 

where  the  subscripts  "B"  and  "S"  indicates  bending  and  the 
transverse  shear  respectively. 

The  bending  and  transverse  shear  stiffness  niatrices  are 
given  as 

[;c/]=//{S^’’B{B^c|ydx  (9) 

0  0 

mb 

Ik*]  =  f f{B^'^G{B^ctydx  (10) 

0  0 


where  E  and  G  are  the  elastic  and  shear  modulii  of  the  beam 
and  the  vectors  {Bg}  and  {Bj}  are  derived  below. 

The  strain-displacement  relationship  for  the  normal  strain 
in  the  axial  direction  is 
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and  the  relationship  for  the  shear  strain  is 


Equations  (11)  2uid  (12)  yield 


dH.  dN,  dK  8fL  8i\Ll 


P  dy  dy  dx  ‘*2  3^  9xj 

The  bending  stiffness  matrix  can  be  obtained  by  carrying 
out  the  integration  in  Eq.  (9)  which  will  result  in 
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For  Eq.  (10)  the  reduced  integration  technique  is  used 
along  the  x-axis  to  prevent  shear  locking  which  occurs  when 
the  ratio  of  beam  length  to  beam  thickness  is  large .  This 
integration  yields  the  transverse  shear  stiffness  matrix  of 
the  form 


[/c/]  = 


Ga 

Ga 

G 

Ga 

Ga 

4b 

4b 

2 

4b 

'  4b 
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The  element  stiffness  matrix,  which  is  obtained  by  adding 
the  bending  and  trouisverse  stiffness  matrices,  can  be 
expressed  in  the  following  form 


[*•]  = 


-a^+aj  a^ 

83-283 

-83-83  -8^ 

ai+2a3  -a^ 

-83-83 

83-83  84 

-^4 

^2 

*4 

-34 

-^2 

-3^-83  a^ 

83+283 

-83+83  -8, 

83-283  -84 

-83+83 

Si +83  84 

^4 

-^4 

^4 

^2 

where 
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The  mass  matrix  can  be  derived  consistently  in  a  similar 
way.  In  this  study  a  diagonal  (lumped)  mass  matrix  was  used 
since  it  takes  less  coirputing  time.  The  mass  matrix  is 


[in*]  =  paJW 


1 

4 

0 

0 

0 

0 

0 


0  0  0  0 

0  0  0 

0  0  0 
2 

0  0  4  0 

4 

0  0  0  4 

4 

0  0  0  0 


0 

0 

0 

0 

0 

1 

2 


(18) 


where  p  is  the  mass  density,  a  and  b  are  the  length  and  height 
of  the  element.  The  element  is  assumed  to  have  a  xmit  depth. 


B.  lAMINXTED  COMPOSITE  BEAMS 

In  the  most  generalized  form  Hooke's  law  can  be  written  as 

using  tensor  notations  [Ref.  4].  Here  represents  the 

elastic  stiffness.  Using  the  symmetry  of  stress  cind  strain 
tensors,  i.e.  and  and  using  a  shorthand  notation 

we  can  write  the  normal  stresses  and  strains  as 
01=011 
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*^2~^22 


‘^2~*^22 

<^3  =  033  EjSPjj 

while  the  shear  stresses  and  strains  are  written  as 


'^4=<^23=<^32 

^4~t23=t32 

t5~t3i5,C33 

^6~t32=t21 

For  general  orthotropic  materials  Hooke's  law  can  be  written, 
in  the  inverted  form,  as 


’ 

^11  ^12  *^13  0  0  0 

/  \ 

S22  S23  0  0  0 

"3 

0 

0 

0 

^4 

0  0 

^4 

^5 

^6 

^5  0 

^5 

^6 

^66. 

(20) 


where  [Sb„]  is  the  conpliance  matrix,  which  is  the  inverse  of 
the  matrix  [C„]  .  For  the  state  of  plain  stress,  Eq.  (20)  is 
reduced  to 


(21) 


Stress-strain  relation  Ceui  be  obtained  by  inverting  the 
relation  given  in  Eq.  (21) 


(22) 
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where  the  reduced  stiffness  terms  are 

l-VjVj 

^22"  1-VjVj  (23) 

n  .  ^2^ 

^  l-VjV^  l-VjVj 

C?66=G6 

In  this  study,  the  width  of  the  beam  is  neglected,  i.e.  it  is 
assvimed  to  have  a  unit  width.  Therefore,  only  the  elastic 
modulus  along  the  beam  axis  is  considered.  Two  different 
laminae  with  four  and  eight  layers  are  used  for  calculations, 
i.e.  [0/90],  and  [0/90/0/90],  where  the  subscript  "s" 
represents  symmetry  with  respect  to  the  middle  axis  of  the 
beam. 

There  are  several  methods  used  in  modeling  laminated 
composite  beams.  In  this  study,  two  different  techniques  are 
used.  The  first  technique  discretizes  respective  layers  in  the 
finite  element  analysis.  As  a  result,  the  total  number  of 
elements  is  proportional  to  the  nximber  of  layers  in  a 
laminated  beam.  Of  course,  static  condensation  can  be 
performed  to  reduce  the  number  of  total  degrees  of  freedom  by 
eliminating  internal  layers  degrees  of  freedom.  This  modeling 
technique  is  computationally  expensive  but  it  can  describe  a 
general  shape  of  deformation  through  the  beam  thickness.  The 
second  technique  uses  one  beam  element  through  the  beam 
thickness  regardless  of  the  niimber  of  layers.  This  technique 
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is  computationally  efficient  but  it  assumes  a  linear 
deformation  through  the  beam  thickness.  The  development  of 
this  technique  is  described  in  the  following  paragraphs. 


Figure  2  :  Laminated  Beam  Element 


Let  Uk  and  v,^  represent  the  axial  and  lateral  displacements 
of  the  kth  layer  and  let  h  be  the  beam  thickness  while  h*  and 
hk*i  r^resent  heights  of  the  top  and  bottom  sides  of  the  layer 
measured  from  the  bottom  of  the  beam  (see  Figure  2)  .  The 
relationship  between  the  layer  displacements  and  the  global 
beam  displacements  can  be  written  as 
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V*.i  =  Vi 

“*%  =  CiUi*+C3U,‘^ 

V*.  ,  =  V, 


(24) 


where  the  constants,  c's,  are 


hr\ 

h 


Cj= 


/,  -  ^**1 


(25) 


This  relationship  can  be  e:q)ressed  in  the  matrix  form  as 


K*1 

n,i 

0  0 
0  0 


0 
0 
0 
0  0 


1  0 


0 

0 

0 


0 

0  0 


K1 

^^2 

(26) 


or 

{d^=[T]{d)  (27) 

where  (d*'}  and  {d}  are  displacement  vectors  of  the  hth  layer 
emd  the  beam,  and  [T]  is  the  transformation  matrix.  Now, 
stiffness  and  mass  matrices  of  a  laminated  beam  element  can  be 
expressed  as 
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(28) 


[T] 

[in*]=j^  [Tl’’[in*]  [T] 

Here,  [k*']  and  [m'‘]  are  the  stiffness  and  mass  matrices  of  the 
kth  layer,  respectively.  Using  this  technique,  a  single  beam 
element  can  include  all  the  layers  of  a  laminated  beam. 


C.  DYNAMIC  ANALYSIS 

In  dynamic  analyses,  displacements,  velocities, 
accelerations,  stresses  eind  strains  need  to  be  calculated  at 
each  time  interval  since  they  all  vary  with  time.  This 
requires  a  time  integration  scheme.  In  this  study,  a  form  of 
central  difference  method  is  used  in  the  dynamic  cuialysis  of 
undanqped  systems. 

The  general  equation  of  motion  at  time  t  can  be  written  as 

(29) 

The  central  difference  method  can  be  applied  in  three  steps. 

A)  Confute  [M]  and  [K] 

B)  Initial  calculations 

1.  WH) 

2.  {u)At  =  {u}  at+{UlQ(At) 

T  'TT  "  (30) 

3.  at  (At) 
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C)  For  each  time  step 


1.  [iC]{u}^) 

2.  {ii).  at={aL  ae  +{llK(At) 

f-j-  t— 5.  t 

2 


(31) 


This  method  is  conditionally  stable  when  the  time  step 
size  At  is  less  than  or  equal  to  the  critical  step  size  [Ref. 
5,6] 


At:SAt  = 

CSt 


(32) 


where  is  the  largest  eigenvalue  of  the  matrix  equation 
(19)  with  no  forcing  function.  When  At>Ater  it  results  in  an 
unstable  nximerical  solution. 


D.  CONTACT- IMPACT  PROBLEMS 

In  the  case  of  partially  delaminated  beams,  as  shown  in 
Figure  3,  the  upper  and  lower  portions  of  the  delaminated 
section  may  overlap  each  other  if  there  is  no  constraint 
applied  on  the  portion  during  the  finite  element  calculations . 
This  fact  applies  to  both  static  and  dynamic  problems .  To 
avoid  this  overlap  the  contact- impact  conditions  are  applied. 


Figure  3  :  Simply  Supported  Beam  with  a  Delaminated  Section. 

The  contact-impact  algorithm  developed  in  [Ref.  7]  is 
applied  to  the  present  beam  element  as  explained  in  the 
following  paragraphs. 

Let  the  letters  i  and  j  represent  the  facing  nodal  points 
in  the  upper  and  lower  portions  of  a  delaminated  section,  and 
let  Vj  euid  Vj  represent  the  vertical  displacements  of  these 
points.  The  procedure  to  apply  the  contact-inpact  condition  is 
explained  below. 

A)  Initially,  the  central  difference  method  is  used  to 
calculate  displacements,  velocities  euid  accelerations  until 
euxy  two  points  contact  each  other  (v^-v^>0). 

B)  VOien  the  two  nodes  penetrate  each  other,  i.e.  v^-v^iO,  the 


16 


release-to-contact  condition  is  applied.  And  the  corrected 
values  of  the  velocity,  acceleration  and  elastic  force  are 
calculated  using  the  following  formula; 


ti  = - li - li 


dj 


T  =T_-- 


(33) 


C)  When  they  are  released  from  each  other,  i.e.  v^-Vj.>0,  the 
contact-to-release  condition  is  applied  and  the  corrected 
values  of  velocity  auid  acceleration  are  computed  using  the 
following  formula; 


2ie 


(34) 


where  the  value  of  a  is  1  for  the  upper  portion  and  2  for  the 
lower  portion  of  the  delaminated  section.  Here  the  subscript 
{  +  )  is  used  for  the  updated  values  while  the  sxobscript  (-1) 
indicates  the  values  taken  at  the  end  of  the  previous  time 
step  and  (-)  indicates  the  values  taken  at  the  end  of  the  last 
iteration  of  the  present  time  step. 

For  static  problems,  the  contact  condition  is  applied 
sinply  by  setting  the  vertical  displacement  of  the  upper 
portion  equal  to  the  vertical  displacement  of  the  lower 
portion  where  they  overlap  each  other. 
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III.  RESULTS  AND  DISCUSSION 


A.  VERIFICATION  EXAMPLES 

The  present  finite  element  model  and  the  computer  code  are 
verified  by  analyzing  three  problems  and  comparing  the  results 
with  exact  solutions.  The  verification  examples  are  static 
deflections  of  both  isotropic  auid  laminated  composite  beams 
and  the  free  vibration  of  beams.  Two  different  boundary 
conditions  are  applied.  One  of  then  is  a  simply  supported  beam 
and  the  other  is  a  cantilever  beam.  When  loading  is  required, 


Figure  4  :  a)  Simply  Supported  Beam  with  Center  Load  b)  Cantilever  Beam  with  End 
Load. 


the  force  is  applied  at  the  center  of  the  singly  supported 
beam  cuid  at  the  end  of  the  cantilever  beam  as  seen  in  Figure 
4. 

Material  properties  of  the  beams  are  as  follows;  elastic 
modulus  E=1,000,000  Pa.,  shear  modulus  G=450,000  Pa.,  density 
p=10  g/cm^,  length  of  the  beam  L=60  m. ,  height  h=0.5  m.  and 
applied  force  P=1  N. 

First  a  free  vibration  analysis  is  performed,  for  both 
simply  supported  and  cantilever  beams.  The  analytical 
expression  for  the  natural  frequencies  is  as  follows 


Yl" 


(35) 


where  y  is  the  mass  per  unit  length  of  the  beam  cund  k^l  are 
defined  values  for  some  beam  types. 

The  results  of  first  three  natural  frequencies  are 
presented  in  Tables  1  euid  2.  As  seen  from  the  results,  the 
finite  element  solution  agrees  with  the  exact  solution  very 
well.  The  mode  shapes  corresponding  to  first  four  natural 
frequencies  are  also  plotted  in  Figures  5  and  6 .  These  are 
ejected  mode  shapes  of  the  simply  supported  and  the 
cantilever  beams  [Ref.  8]. 
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TABLE  1  :  NATURAL  FREQUENCIES  OF  SIMPLY  SUPPORTED  BEAM 


'''V"::;;',; 

Esedct  Solution 

0.00395 

0.0158 

0.0356 

P8W  solution 

0.00395 

0.0158 

0.0356 

- 

- 

- 

TABLE  2  :  NATURAL  FREQUENCIES  OF  CANTILEVER  BEAM 


1 

% 

Q3 

1  Exact  solution 

0.00141 

0.00883 

0.0247 

0.00141 

0.00883 

0.0247 

Hi 

- 

- 

- 

Then  static  defl^'ctions  of  beams  are  computed.  The 
analytical  solutions  for  isotropic  beams  are  [Ref.  9] 

where  Eq.  (36)  is  for  a  simply  supported  beam  and  Eq.  (37)  is 
for  a  cantilever  beam,  here  "I"  refers  to  the  moment  of 
inertia  and  L  is  the  length  of  the  beam.  The  results  of  both 
exact  solution  and  finite  element  solution  are  plotted  in 
Figure  7 . 
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mode  1 


mode  2 


Figure  5  ;  Mode  Shapes  of  the  Simply  Supported  Beam 
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mode  1  mode  2 


Figure  6  :  Mode  Shapes  of  the  Cantilever  Beam 
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Figure  7  ;  Static  Deflections  of  Isotropic  Beams;  a)  Simply  Supported  Beam , 
b)Cantilever  Beam 
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Finally  a  static  cuialysis  of  a  laminated  composite  becun 
with  simple  supports  is  performed.  The  analytical  expression 
for  the  deflection  of  a  simply  supported  beam,  which  is 
subject  to  three  point  bending,  as  shown  in  Figure  8,  is  as 
follows  [Ref.  10] 

PL^x 

- TT 

48B*J 

cind  the  equivalent  bending  stiffness  is  defined  as 

(39) 

Jlc=l 

here  is  the  effective  bending  modulus  of  the  beam,  Ei’'  is 
the  modulus  of  the  Icth  layer  relative  to  the  beam  axis,  l’‘  is 
the  moment  of  inertia  of  the  Jcth  layer  relative  to  the 
mic^lane,  and  N  is  the  number  of  layers  in  the  laminate.  The 
vertical  displacements  obtained  from  the  exact  solution  and 
the  finite  element  solution  are  plotted  together,  for  four  and 
eight  layered  beams  and  for  two  different  ratios  of  elastic 
modulii,  in  Figures  9  and  10. 
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Figure  8  :  Laminated  Beam  with  Three  Point  Bending  (Simply  Supported). 

In  the  static  problems,  as  seen  from  the  figure  the 
exact  solutions  and  the  finite  element  solutions  agree  very 
well.  From  these  verification  exaztples,  it  can  be  shown  that 
the  beam  element  provides  accurate  results  compared  to  exact 


solutions . 


displacement,  m. 


a) 


length,  m. 


Figure  9  ;  Static  Deflections  of  the  Laminated  Composite  Beam  with  4  Layers 
a)forEo/E^=20,  b)  for  Eo/Ejo^lOO 
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cNsplaceinenl,  m. 


a)forE,/EM«20,  b)  for  Eq/Em»100 
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B.  NOMBRZCAL  RESULTS 

Dyncunic  analyses  are  performed  for  a  siirply  supported 
composite  beam  with  four  layers,  i.e.  [0/90],,  and  with 
different  crack  sizes  using  contact -impact  conditions. 
Calculations  without  contact-impact  conditions  are  also  made 
and  compared  to  the  case  computed  with  contact-iirpact 
conditions.  An  impulsive  load  of  a  very  short  duration  is 
applied  at  the  center  of  the  upper  portion  of  a  delaminated 
section  and  the  displacement  response  of  the  beam  is  recorded 
for  a  period  of  time. 

First  displacements  at  the  mid-points  of  upper  and  lower 
portions  of  a  large  crack  are  calculated  and  plotted  in  Figure 
11  without  using  contact -inpact  conditions,  where  crack  length 
is  one  third  of  the  beam  length.  As  seen  from  the  figure, 
these  two  points  overlap  each  other  which  is  not  admissible. 
Then  for  the  same  crack  size,  contact -impact  conditions  are 
applied  and  displacements  of  the  beam  is  calculated  and 
plotted.  Figure  12  shows  the  displacements  of  the  beam  just 
before  contact  while  the  displacements  at  the  time  of  contact 
is  shown  in  Figure  13 .  The  displacements  just  after  contact  is 
shown  in  Figure  14.  Also  the  responses  of  the  mid-points  and 
quarter -points  are  plotted  in  Figxires  15  and  16,  respectively, 
where  gueurter -points  represent  the  points  which  are  located  in 
the  middle  between  the  crack  tip  and  the  crack  mic^oint  for 


both  upper  and  lower  portions.  Later,  the  crack  size  is 
reduced  to  one  thirtieth  of  the  beam  length  and  the 
displacements  of  mid-points  and  quarter-points  are  plotted  in 
Figures  17  and  18,  respectively.  As  seen  in  these  results,  the 
upper  and  lower  portions  never  penetrate  each  other .  This 
indicates  that  the  application  of  contact- impact  conditions 
works  well  in  simulating  delaminated  beams.  Finally,  the 
response  of  a  beam  without  a  crack  is  computed  and  plotted  in 
Figure  19  for  comparison  purpose.  Responses  of  all  these  cases 
are  different  each  other.  They  may  give  cui  idea  how  cui 
exnbedded  crack  affects  the  behavior  of  beams. 
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Figure  11  ;  Displacements  of  Mid-Points  without  Contact-Impact  Conditions 
(Large  Crack  Size  ) 
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Figure  12  ;  Displacement  of  Beam  (before  Contact) 


Oisplaceroent,  m. 


t-2.923  sec. 


Length,  m. 


Figure  14  :  Displacement  of  Beam  (after  Contact) 
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Figure  15  ;  Displacements  of  Mid«Points  with  Contact-Impact  Conditions 
(Large  Crack  Size) 
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displacement,  m. 


Figure  17  ;  Displacements  of  Mid-Points  with  Contact-Impact  Conditions 
(Small  Crack  Size) 
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Figure  18  ;  Displacements  of  Quarter-Points  with  Contact-Impact  Conditions 
(Sn^i  Crack  Size) 
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IV.  CONCLUSIONS  AND  RECOMMENSATIONS 


The  beam  element  used  in  this  study  is  useful  to  model  a 
crack  w’  '  contact-impact  conditions.  Application  of  contact- 
impact  -iditions  constrains  the  vertical  displacements  at  the 
crack  in  a  way  that  the  points  on  lower  and  upper  portions  of 
the  crack  hit  each  other  without  penetration, 

Conputational  time  depends  on  the  nvimber  of  degrees  of 
freedom.  The  technique,  which  models  the  laminated  composite 
beam  using  a  single  beam  element,  reduces  the  degrees  of 
freedom  to  save  time  during  calculations. 

When  an  impulsive  load  of  a  short  duration  is  applied,  the 
beam  responds  in  different  ways  depending  on  whether  there  is 
a  crack  inside  it  and  what  the  crack  size  is.  Examining  these 
results  gives  an  idea  about  existence  of  an  embedded  crack  as 
well  its  size  and  location. 

For  further  studies,  delaminated  sections  may  be  refined 
to  see  the  minimum  crack  size  that  can  be  detected  using  this 
procedure.  Also,  the  Fast  Fourier  Treuisform  can  be  applied  to 
measure  the  frequency  response  of  a  delaminated  beam. 
Comparing  the  natural  frequencies  may  give  an  idea  about  the 
crack  size.  In  this  stu<^,  a  crack  is  assumed  to  be  along  the 
middle  axis  of  the  beam  and  the  load  is  applied  at  the  center 
of  the  beam.  Locations  of  cracks  and  loads  may  be  changed  and 
the  responses  can  be  calculated  for  different  cases. 
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